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ABSTRACT 

Stellar mass compact objects in short period (P < 10 3 s) orbits about a 10 4 5 -10 7 5 M© massive 
black hole (MBH) are thought to be a significant continuous-wave source of gravitational radiation 
for the ESA/NASA Laser Interferometer Space Antenna (LISA) gravitational wave detector. These 
extreme mass-ratio inspiral sources began in long-period, nearly parabolic orbits that have multiple 
close encounters with the MBH. The gravitational radiation emitted during the close encounters may 
be detectable by LISA as a gravitational wave burst if the characteristic passage timescale is less than 
10 5 s. Scaling a static, spherical model to the size and mass of the Milky Way bulge we estimate an 
event rate of ~ 15 yr _1 for such burst signals, detectable by LISA with signal-to-noise ratio greater 
than five, originating in our Galaxy. When extended to include Virgo Cluster galaxies our estimate 
increases to a gravitational wave burst rate of ~18 yr _1 . We conclude that these extreme mass-ratio 
burst sources may be a steady and significant source of gravitational radiation in the LISA data 
streams. 

Subject headings: black hole physics — Galaxy: nucleus — gravitational waves — stellar dynamics 



1. INTRODUCTION 

The inspiral of compact objects onto massive 
black holes (MBHs) in galactic nuclei is an antic- 
ipated important gravitational radiation source for 
the ES A/NASA Laser Inter f erometer Space Antenna 
(USA) (iHils fc Benderl IT9951 iSieurdsson fc Reesl Il997t 
GlamDcdakis & Kenncfick 2002; GlamDcdakis ekjiTJ 
20021 llvanovl 120021 iFreitad l2003 lOair et all | 2004f 
Danz maim fe Rudieerl l2003t jSumner fe Shaul 120041 
Jennrichl 12005: ISetoetall l200l|) " These extreme mass- 



ratio inspiral (EMRI) sources will complete > 10 4 orbits 
during the LISA mission lifetime (~ 5 yrs) and it is the 
accumulation of the signal power, emitted continuously 
with frequency > 10 -3 Hz, that makes them visible to 
LISA. 

However, before emerging as EMRIs these objects were 
on long-period, nearly radial orbits that were "captured" 
by the MBH after the orbital energy was reduced by a 
series of short, intense bursts of gravitational radiation 
emitted during each close encounter with the MBH. If 
the encounter timescale is less than ~ 10 5 s the gravita- 
tional wave burst will be in the LISA band and, if strong 
enough, detectable by LISA. Here we describe a prelim- 
inary analysis of this extreme mass ratio burst (EMRB) 
phenomenon, characterizing the detectability and esti- 
mating the rate of EMRB events in LISA. Our estimates 
show that the Milky Way may be responsible for ~ 15 
observable EMRB events per year, and that the Virgo 
Cluster galaxies may be responsible for an additional ~ 3 
events per year, making EMRBs a significant source of 
gravitational waves for LISA. 

2. INSPIRALS VS. BURSTS 

To illustrate the evolution from EMRB to EMRI con- 
sider a 0.6 M© white dwarf in a nearly radial Keplerian 
orbit about a MBH of mass M. ~ 10 6 M . If the ini- 
tial orbit has an apocenter of 100 pc and a pericenter 
of 0.2 AU (corresponding to 10 Schwarzschild radii for 



the MBH) , then during the first pericenter pass the orbit 
loses enough energy via gravitational radiation that the 
next apocenter is ^60 pc, while the change in pericenter 
is negligible. With each subsequent pericenter pass, a 
strong burst of gravitational radiation reduces the sub- 
sequent apocenter and the orbit becomes more circular. 
LISA is sensitive to gravitational waves in the approxi- 
mately 10~ 5 to 10 -1 Hz band, so we classify as EMRBs 
those systems with orbital periods greater than 10 5 s 
and with pericenter passage timescales of less than 10 5 s. 
Conversely, EMRIs have orbital periods less than or on 
order 10 3 s and radiate continuously in the LISA band. 

EMRIs and EMRBs sample different components of a 
galactic nucleus. Every EMRI source was, in its past, a 
possible EMRB source: i.e., it was on a highly eccentric 
orbit that, owing to gravitational wave emission, circu- 
larized and decayed until it entered the LISA band as 
a continuous source. However, if the MBH encounter 
timescale is too long, the orbit may become an EMRI 
without any significant radiation burst in the LISA band. 

Not all EMRBs evolve to become EMRIs, either. An 
object on a highly eccentric orbit may be scattered by 
other stars after encountering the MBH, or may plunge 
directly into the MBH. In addition, while massive main 
sequence stars will be disrupted long before appearing as 
EMRIs, the pcriapsis of an EMRB orbit can be much 
greater than the tidal disruption radius for low-mass 
main sequence stars. For example, an M2V star could 
pass within 0.65 AU of the Milky Way's MBH without 
disruption and, if on a nearly radial orbit so that its 
pericenter velocity was highly relativistic, may radiate a 
gravitational wave burst with a characteristic frequency 
of several mHz. 

3. STELLAR MODEL 

To evaluate the EMRB event rate we begin with a 
model of our Galaxy center. The Milky Way bulge 
consi sts of a M. = 3.7 x 10 6 M© MB H llOhez et al 



12005]) embedded in a stellar ellipsoid (e.g. iBinnev et al 



119971 IStanek et all 119971 ILopez-Corredoira et all 12000(1 

of mass (1.3-2) x 10 10 M^ le.g. lGerhardll2002tt. with axis 
lengt hs of ~ 1.8 : 0.7 : 0.5 kpc ijBissantz fc GerhardI 
2002), and with a cuspy density profile that tend s to- 
ward p ~ r -1 - 8 as r -» l|Matsumoto et al.lll982fl . Al- 
though the Milky Way is a barred galaxy, we treat the 
bulge a s spherical with a dens ity profile described by an r) 
model llTremaine et al.l fl994). Since rj models are static, 
spherical, and isotropic their distribution function de- 
pends only on energy, making the model an ideal first 
choice to analytically estimate the number of bursts in 
the bulge. We choose ij = 1.25 to match the observed 
slope of the Milky Way inner cusp and set the effective 
radius equal to 2 kpc, total stellar mass to 2 x 10 10 M , 
and MBH mass to 4 x 10 6 M Q . 

Embedding an r\ model with a pre-existing central 
MBH leaves the stellar density profile unchanged, but 
changes the cluster potential so that $.(r) = $*(r) — 
M,/r, where <£. is the potential of the MBH-embedded 
model and is the stellar potential of an n model. The 
new potential naturally adjusts the distribution function, 
which can be obtained from a sph erical, isotropic den- 
sity p rofile via Eddington's formula ( Binn ev fc Tremainel 
Il987f) . which for a MBH-embedded model can be solved 
analytically near the MBH: 
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where e is the absolute value of the energy per unit 
mass. When normalized to the total mass of the model, 
f(e)drdv is the mass contained in the local volume ele- 
ment centered on r, v. 

4. BURST WAVEFORMS 

Orbits that may generate observable EMRBs are char- 
acterized by large orbital periods (~ 1 — 10 9 yrs) and 
high eccentricity. Objects on these orbits may be highly 
relativistic during periapsis passage, with velocities an 
appreciable fraction of the speed of light (v p ~ 0.3c). 
The gravitational wave emission from such a "gravita- 
tional bremsstrahlung" event emerges as a burst, lasting 
as long as the MBH encounter and beamed in the direc- 
tion of the secondary's velocity at periapsis. For this pre- 
liminary exploration we treat trajectories near periapsis 
as Keplerian, 1 with the same v p , r p , and mass enclosed 
within r p as dictated by the rj model, and we calculate 
the gravitational radiation as if it were well described 
by the quadrupole formula. A more detailed treatment 
would include the effect of beaming, the radiative contri- 
butions from terms in v/c higher than quadrupole, and 
take into account the actual relativistic trajectory, all 
of which may be impo rtant for traject ories with small 
r p /M. or large v/c Ccf. lGair et aDl2005^ . 

From the (quadrupole) waveforms we calculate the 
signal-to-noise ratio (SNR) in the LISA detector assum- 
ing the source is at a distance of 8 kpc. In general, the 
square of the SNR is given by, 
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1 Although we only considered elliptical systems within the rj 
model, it is worth noting that the Keplerian equivalent orbit near 
the MBH is not necessarily elliptical. 



where h(f) is the Fourier transform of the gravitational 
wave signal projected onto the LISA detector and S n (f) 
is the one-sided instrument noise power spectral density. 
For our purposes, the values of S n (f) ar e taken from 
the O nline Sensitivity Curve Generator ijLarson et alJ 
2000) with the standard LISA settings and the inclu- 
sion of a white dwar f background "noise" contribution 
lIBender fc Hilslll997h . 

To gain insight into the magnitude and scaling of the 
SNR, equation [21 can be approximated by taking ad- 
vantage of the simple burst structure in the waveforms, 
which allows us to quickly approximate the integral and 
Fourier transform, 
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where R is the distance to the EMRB, M* is the mass 
of the secondary, and /* = v p /r p is the inverse of the 
burst duration. The value of /* can be viewed as a 
characteristic frequency for an EMRB event, but these 
events are extremely broadband in nature. Note that the 
noise spectral density scales as f~ 4 for frequencies below 
1 mHz and the pre-factor accounts for the white dwarf 
binary background. 

5. EVENT RATES 

To arrive at an EMRB event rate estimate, we divided 
our r\ model into five million elements in (r, v) phase 
space. We chose the element centers uniformly spaced 
in log(r) and log(u) with r G [7.6 x 10~ 7 , 2 x 10 4 ] pc and 
v € (0, 2 x 10 5 ) km s _1 . For each phase-space element we 
calculated /(e), the average orbital period P or b, periapsis 
r p , and velocity at periapsis v p . 

To arrive at the number of compact objects and low 
mass main sequence stars in our model we do the fol- 
lowing. Assuming all of the mass density in our model 
results from a single burst of star formation 10 Gyr ago, 
and that the number of sta rs per unit mass follows a 
Kroupa IMF 2 ijKroupal 200lj), we converted /(e) into the 
total number of stars per mass per phase space element. 
This number can be less than one per phase space ele- 
ment. To determine the number of main sequence stars 
per phase space element we calculated the maximum stel- 
lar mass that would survive tidal disruption by the MBH 
(|Sridhar fc Tremainel[i992l) and integrate the IMF up to 
this mass, setting the main sequence stellar mass equal to 
the average mass that survives tidal disruption. We de- 
termined the white dwarf population by assuming white 
dwarfs formed from all stars between 1 — 8 M and the 
final white dwarf mass is 0.6 Mq. Similarly, neutron 
stars were assigned masses of 1.4 M Q and form from 
stars 8 — 20 M Q , while black holes are 10 M objects 
formed from 20— 120 M stars. Our calculation neglects 
the effects of mass segregation, multiple episodes of star 
formation, and forces all stars between 1 — 120 M Q to be 
a stellar remnant. 

Only a small fraction of the initial phase space con- 
tain EMRBs. To isolate the EMRB orbits we made 

2 We varied the IMF choice to include Salpeter and Scalo IMFs; 
this did not change our rate significantly. 
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Fig. 1. — Orbital parameters of potential EMRB sources. Here e 
is the eccentricity of the orbit within the total bulge potential, 
and R p is in units of Schwarzschild radii. Typical sources are 
more eccentric and have larger pericenter distances than the av- 
erage EMRI. Overlaid is the subset of white dwarf orbits that are 
observable with LISA at a SNR > 5. 



a series of cuts. We immediately discarded those or- 
bits that were either unbound or had periapsis within 
four Schwarzschild radii. We excluded orbits with P or b 
greater than the relaxation time, t re i ax = v 2 / D(Av 2 ), 

where D(Av 2 ) is the Fokker-Planck diffusion coefficient 

ijBinnev k, Tremaindll987l eq. 8-68). We also excluded 
orbits with P OT b less than 1 5 s, which would yield a con- 
tinuous signal in the LISA band and are more accurately 
classified as EMRIs. We also excluded orbits that radiate 
so strongly that they plunge into the MBH faster than a 
dynamical time, where the plunge timescale T p i unge is 
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Here r p ^ is the initial pericenter distance and ej is the 
initial eccentricity. Finally, we excluded all orbits with an 
encounter timescale At — r p /v p > 10 5 s, as these passes 
emit radiation outside the LISA band. The remaining 
stars constituted our potential EMRB population. 

Figure ^ shows the phase space remaining after these 
cuts; ~ 1, 800 phase space elements remained, which does 
not imply 1800 distinct objects; in fact there are typically 
10~ 7 white dwarfs per drdv in this region. Overlaid in 
the same figure is the subset of the white dwarf phase 
space with SNR > 5. 

The total event rate v is the sum over all orbits with 
LISA SNR greater than 5: 
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For our Milky Way bulge model, we find javd = 
z^ns = 0.1 yr -1 , v B n = 0.06 yr -1 , and ^ L MMS,survive 
1 2 yr -1 . These event rates imply that EMRBs may be 
an important, heretofore unrecognized source of gravita- 
tional waves in the LISA band. The left panel of figure |5] 
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Fig. 2. — Left: Accumulative event rate as a function of the 
orbital period. MS = main sequence; WD = white dwarf; NS 
= neutron star; BH = black hole. Right: Number of black hole 
sources versus signal to noise ratio. All sources above the light grey 
section have SNR > 5. Note the small fraction of EMRBs with 
SNR> 10 4 (dark grey), which would be visible to Virgo distances. 



shows the accumulative event rate as a function of orbital 
period. Most EMRBs originate from orbits that would 
produce multiple bursts of radiation per year in the LISA 
band. The right panel of figure |2 shows the number of 
black hole EMRBs as a function of SNR; note that some 
have SNRs high enough to be seen out to Virgo distances. 
If we repeat our analysis with the EMRB distances set 
to 16 Mpc and multiply by the number of galaxies in the 
cluster with suitable mass MBHs, we find a Virgo event 
rate of ~3 yr -1 , all due to encounters of low mass black 
holes with central MBHs. 

6. DISCUSSION 

EMRBs may comprise a significant new low frequency 
gravitational wave source, rivaling the total number of 
EMRI event s visible in an y year of LISA observation. 
For example. iFreitad l|2001|) estimates that O(0. 1) white 
dwarfs per year are in the EMRI phase in the Milky 
Way. We find 0(1) EMRBs per year from close white 
dwarf encounters. The relaxation requirements on EM- 
RBs are much less strict: EMRIs require thousands of 
undisturbed orbits to accumulate a SNR large enough 
to be detected (T p [ unge > (1 — e)T re i ax ), while EMRBs 
merely need to pass by the MBH before scattering via 
2-body relaxation. 

Our estimated Milky Way EMRB event rate of ~ 
15 yr -1 may appear high in light of other processes, such 
as the white dwar f capture rate , which estimates place 
0(1O~ 7 ) per year l)Freitadl2001j) ; however, the greatest 
contribution to the event rate we calculate arises from 
sources that burst multiple times a year, yet individu- 
ally have a small phase space density. In other words, 
our rate does not imply that 15 Mq of material interact 
with the black hole per year; rather, given these EMRB 
orbits we calculate that Milky Way black hole grows by 
less than 10 4 M Q over a Hubble time via EMRB decay, 
with many EMRBs bursting hundreds of times before 
falling into the MBH. Our simple stellar model passes 
ot her sani t y-chec ks as well: given the capture definition 
of IFreitad" (|2001fl . we expect a white dwarf capture rate 
of 0(1O -6 ) per year and find, in our calculations, a tidal 
disruption rate of 0(1O -5 ) p er year, which is con sistent 
with observational estimates ijDonlev et alJ l2002'). 

We have introduced EMRBs as the evolutionary pre- 
cursor to EMRIs; many EMRBs will inspiral and even- 
tually evolve into EMRIs as gravitational radiation car- 
ries away energy and angular momentum. This implies 
that the lifespan of the typical EMRB phase may only 
be 0(1O 7 ) — 0(1O 9 ) years after galaxy formation or a 
star formation episode. We have assumed, however, that 
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EMRB phase space is always occupied. For there to 
be a significant population of EMRBs in the present- 
day galaxy sample, there must be a refilling mecha- 
nism. There are several mechanisms that may refill 
the burst reservoir. One such mechanism, the dynam- 
ical migration of a stellar cluster, may have recently 
left its mark in the inner 0. 04 pc of our gala xy in the 
form of the young S stars llKim et al.ll2004 |). Other 
mechanisms, such as triaxiality lIMerritt fc Poor] 120041 
iHollev-B ockelmann fe Sigurdsso nl 120061) or resonant re- 
laxation n|Ho^nilm*7r^kx£uidern'2006f) . can act to refill 
the loss cone on timescales much shorter than the two- 
body relaxation time. 

Since we were primarily concerned with determin- 
ing whether EMRBs were an overlooked class of LISA 
sources, we considered the simplest possible model for the 
Milky Way bulge, the radiation properties, and the abil- 
ity of gravitational wave analysis techniques to detect the 
signal. Under these approximations the predicted event 
rate implies that bright EMRBs are so numerous that the 
local universe will generate a background of gravitational 
wave bursts with a mean rate of once per three weeks. 
Changing the assumptions, however, can alter this rate 
by several orders of magnitude; modeling the bulge as a 
bar, for example, can increase the capture rate by two 
order s of magnitude ( Hollcv-Bockcl mann fc Sigurdssora 
2006), while mass-segregation can decrease the EMRB 
rate from lo w mass main sequence stars by two orders of 
mag nitude l|Freitagll200l . 



The relatively large event rate calculated here suggests 
that EMRBs may in fact be a significant source for the 
LISA detector and indicates that further, more accurate 
modeling needs to be undertaken. To better character- 
ize EMRBs and study how they can constrain models of 
galactic nuclei, future work needs to include a more accu- 
rate model for the gravitational wave emission, a better 
treatment of the orbit in the relativistic regime, more 
realistic time-dependent galaxy and star formation mod- 
els, and a proper treatment of stellar dynamics, including 
mass segregation. In addition, the potential presence of 
these burst gravitational wave sources introduces a new 
set of data analysis challenges, associated with the detec- 
tion and characterization of the EMRB signal from LISA 
data. As daunting as the prospect of these challenges are 
their potential payoff: the signal from EMRBs probe the 
central region of our Galaxy, and galaxies in the Virgo 
Cluster, at scales less than 10 -5 pc, far greater resolution 
than we can ever hope to achieve electromagnetically. 
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